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O ■ Abstract 

The extra contributions due to an extra generation of vector-like quarks to the isospin symmetry 
Q-i, breaking observable associated with the radiative B K*^ decay is obtained. It is shown that 

' this additional contribution is sensitive to the nonunitarity parameter U^^, which is a measure 

■ of the the strength of the non-zero tree-level flavor changing neutral current in this model. The 

. significance of this result is that, once accurate experimental results on the isospin asymmetry 

becomes available, one can constrain the U^^ independent of the mass of the extra quarks and so 
a much more stringent acceptable model parameter space could be obtained. 

PACS numbers: 13.20.He, 12. 39. St, 12.15.Mm 
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I. INTRODUCTION 



The Standard Model (SM) of electroweak interactions had been quite successful in de- 
scribing the experimental data up to recent past. However, there are now some indications 
of possible inconsistencies between this model and experiment. In order to improve our 
understanding of what lies beyond the SM, we need to make progress on two fronts: 1) 
improve the accuracy of the experimental data 2) reduce the uncertainties in the theoretical 
predictions of the physical observables based on the SM. The hope is to find the deviations 
from the SM on as many front as possible so that to be able to rule in or out the various 
scenarios for the new physics. 

One of the main sources of theoretical error in making predictions involving hadrons is 
the nonperturbative nature of the strong interactions. Since the weak interactions involving 
quarks is one of the main venues to test the SM, it is quite essential to be able to deal with 
this issue. Due to confinement, the physical states that appear in the matrix element of 
the quark weak currents are hadrons, which are bound states of quarks, rather than free 
quarks themselves. The form factors that parameterize these matrix elements characterize 
our lack of knowledge of the strong force in the nonperturbative regime and can not be 
evaluated analytically from first principles. However, with the help of the fast computers, 
lattice QCD has been able to make significant progress in numerical calculation of some 
of the nonperturbative parameters. Nonetheless, to reduce the theoretical uncertainly, one 
solution for the time being is to look for observables which are not too sensitive to form 
factors and other nonperturbative parameters. 

The other criteria for an interesting process for testing the SM is the sensitivity to new 
physics. Usually, one looks for a process which receives vanishingly small contribution from 
the SM and therefore, any significant deviation from zero in the experimental measurement 
could be attributed to the new physics beyond the SM. On the other hand, experimental 
confirmation of the null result to a good precision should rule out a good number of new 
physics scenarios if the process satisfies the above criteria. 

The radiative inclusive decay B — > Xs'y, which is due to the underlying fiavor changing 
neutral current (FCNC) quark transition b —>■ sj, and its exclusive mode B K*j have 
proven to be quite important processes for examining the SM and new physics since they 
were first observed by CLEO. The most recent experimental data on these transitions are 
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as follows 



HQS: 



Br{B+ K*+j) = < 



(4.01 ± 0.21 ± 0.17) X 10-5 (BELLE) 

(3.92 ± 0.20 ± 0.24) x 10"^ (BABAR) 

(4.55l°;^2 ^ 0.34) x 10"^ (CLEO) 

(4.25 ± 0.31 ± 0.24) x 10"^ (BELLE) 

(3.87 ± 0.28 ± 0.26) x 10^^ (BABAR) 

(3.76^°:| ± 0.28) X 10-5 (CLEO) 



(1) 



isQ 



le current PDG average for the inclusive mode from CLEO and the BELLE measurements 



Br(B X,7) = (3.3 ± 0.4) x W 



(2) 



This radiative transition is dominated by the one-loop quantum effect in SM which in 
effective Hamiltonian language is represented by the magnetic dipole operator as illustrated 
in fig. n Since the exotic particles that appear in various extensions of the SM can appear 
as quantum fluctuations in the loop, this rare B decay has been at the focus of many 
theoretical investigations to constrain different new physics scenarios. The inclusive decay, 
though easier to handle theoretically, is much more difficult to be accurately measured. To 
leading order in as, one can equate the B Xg'y equal to the underlying quark transition 
b ^ s'y (fig. P). The systematic approach to decays at energies much below the VT-boson 
mass, like B-meson decays, is through the effective Hamiltonian. Here the heavy fields, like 
W and top-quark fields, are integrated out of the electroweak lagrangian and the effective 



Y 



FIG. 1: The effective 5 — > 57 transition. 
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Hamiltonian is written in terms of a series of operators with increasing mass dimensions. Of 
course, practically, only the lowest and next-to-lowest mass dimension operators contribute 
to the transitions. These operators are listed as follows: 

• Current-current operators: 

Olc = qa1^LCaCp-i^,Lhp , O2C = QaYLcpCp-ff^Lha . 

• QCD penguin operators: 

O3 = qaYLba q'plt^Lq'p , ^4 = q^l^Lbp ^ q'plt^Lq'^ , (3) 

q' q' 

O5 = qal^Lba ^'p^M ' ^6 = qal^Lbp ^ q'pl^.Rq'a ■ (4) 

q' q' 

• Electroweak penguin operators: 



q'=u,d,s,c,b 
q'=u,d,s,c,b 

3 
3 

Oio = -qal'^Lbp Y eg'g^7/.^9a • 

q'=u,d,s,c,b 

• Electro- and chromo-magnetic operators: 

g 

= -^qaCr^^irribL + msR)baF^,y , (6) 
Og = ^^q.a^\m,L + m,R)T:pbpGl, . (7) 

In the above equations q = s or d and L{R) = are the projection operators. As 

mentioned above, has the main contribution to the b —>■ s'y transition in the leading order. 
At the next-to- leading order in as, besides the electromagnetic penguin operator, the current- 
current operators Oi and O2 (fig. El) and the chromomagnetic penguin operator Og(fig. 
contribute as well. This leads to the SM prediction of Br{B X^'-y) = (3.6 ± 0.30) x 10~^ 
and consequently provides constraints on the new physics beyond the SM.j 



The exclusive B — ^ K*'y decay is more difficult to calculate theoretically due to the 
nonperturbative nature of the hadronic matrix element of the above operators. These matrix 
elements are parameterized in terms of a number of form factors. One could use certain 
approximate symmetries of the QCD, like heavy quark symmetry, to reduce the number of 
these parameters. But in any case, the calculation of these parameters introduces model 
dependence and theoretical uncertainty in our predictions. 

In Ref. jsj, a model independent approach based on QCD factorization in the leading 
power of h^Qco/^B has been proposed. Here, the operators that were mentioned previously 
in connection to the inclusive decay contribute to the factorizable term to the leading (fig. |3} 
and next-to-leading order in (fig- EI)- On the other hand, there are new operators involving 
the light quark in the B meson that contribute to the nonfactorizable term as illustrated 
in fig. ini). Also, there are annihilation diagrams contributing to the nonfactorizable decay 
of the charged B mesons (fig. [7j), however these transitions are CKM as well as power 
{^qcd/^b) suppressed. The nonperturbative QCD enters through a number of universal 
functions, a form factor for B K* transition and light-cone distribution amplitudes for 
B and K* mesons. In fact, due to the sensitivity of the branching ratio to the flavor of the 
light quark through the nonfactorizable terms in the ones depicted in fig. El one obtains a 
slightly different branching ratio for neutral and charged B mesons. Taking the appropriate 
value for the input parameters, authors in Ref. 0] obtain Br{B^ K*^'~f) = 7.45 x 10~^ 
as compared to Br{B^ K*^-f) = 7.09 x 10~^. 




FIG. 2: 0{as) current-current operators contribution to 6 — > ,57. Cross marks are the alternative 
locations for the coupling of the emitted photon. 
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II. ISOSPIN SYMMETRY BREAKING IN S ^ K*-^ 

Besides the branching ratio, one could associate another observable, the isospin symmetry 
breaking, to the radiative B —>■ K*'-f decay. This asymmetry is define as follows: 

^ r(^" ir*°7) - T{B- ^ K*-^) 

with Ao+ obtained from eq.® by using the charge conjugate modes. Ao± could prove to 
be important observables for examining the SM as well as discriminating between various 
new physics scenarios y. The data from Belle Q and Babar[2] point to isospin asymmetries 
of at most a few percent and consistent with zero within the experimental error: 

Ao- = +0.051 ±0.044(stat.) ±0.023(sys.) ± 0.024(i?+/°) (Babar) , (9) 
Ao+ = +0.012 ±0.044(stat.) ±0.026(sys.) {Belle), (10) 

where the last error in eq.Q is due to the uncertainty in the ratio of the branching fractions of 
the neutral and charged B meson production in T(4S') decays. As we noticed in the previous 
section, this asymmetry in the SM is due to the non-spectator contributions and has been 



estimated within the QCD factorization approach in Refs. [y| and 



a 



and the perturbative 



QCD method in Ref. j^. In general, to the lowest order in and the small parameter 
^qcd/itT'B, current-current, QCD penguin and chromomagnetic penguin operators make 
the dominant contributions to the isospin symmetry breaking within the SM. The role of 
the electroweak penguin operators are more or less the same as their QCD counterparts, 
however, since their corresponding Wilson coefficients are significantly smaller, they can 
be safely ignored. The more accurate measurement of the isospin asymmetry in the near 




FIG. 3: 0{as) chromomagnetic penguin operator contribution to 6 — > 57. Cross mark is the 
alternative location for the coupling of the emitted photon. 
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future and a better understanding of the SM prediction for this observable should provide a 
sensitive testing venue for possible models of new physics. One such model is the extension 
of the SM with an extra generation of iso-singlet quarks jl^. 
next section. 



Ill ] which is explained in the 



III. VECTOR-LIKE QUARK MODEL (VQM) 

In this model, the gauge structure of the SM remains intact except for an additional pair 
of iso-singlet quarks, which we denote them by U and D. The difference between these new 
quarks and ordinary quarks of the three SM generations is that, unlike the latter ones, both 
left- and right-handed components of the former quarks are SU{2)l singlets. Therefore, the 
Dirac mass terms of vector-like quarks, i.e. 

muiULUR + URUL)+mDiDLDR + DnDL) , (11) 

are invariant under electroweak gauge symmetry. However, the masses of the ordinary quarks 
arise from their gauge invariant Yukawa couplings to an iso-doublet scalar Higgs field as 
follows: 

-ff't-iM - fS'l-lKi + H£: , (12) 

q q 




FIG. 4: Leading order factorizable contribution to i? — > K*'y. 
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where i,j — 1, 2, 3 is covering the three generations of the regular quarks, and the doublet 
of fermions is defined as 




At the same time, additional SU{2)l invariant Yukawa couplings between vector-like and 
ordinary quarks, in the form 

-fii'lD^'t' - J"'l'lUn4> + H.C. , (14) 




q q 




FIG. 5: Ncxt-to-leading order factorizable contributions to 5 — > K*^. Cross marks are the alter- 
native locations for the coupling of the emitted photon. 
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leads to mixing among 4 up- and down-type quarks of the same charge. As a result, after 
spontaneous electroweak symmetry breaking due to {(p) — v ^ 0, we obtain the following 
mass terms: 

+ ulMfu^^ + H.C. , (15) 

where and being 4x4 mass matrices and a, (3 — 1..4 cover ordinary and vector- 
like quarks. In general, the mass matrices are not diagonal and unitary transformations 
from weak to mass eigenstates are necessary to achieve diagonalization. Denoting the mass 
eigenstates with u'^ ji and d'^ j^, we have 

„,a _ /IM al3 I /3 ,a — Ad "1^ J' /-if^N 

where the unitary transformation matrices A^^j^ are chosen such that Af'^MdAj^ and 
^ are diagonal. The interesting property of the VQM is that the transformations 
(6) lead to inter-generational mixing among quarks not only in the charged current sector 
but also in the neutral current interactions. This is due to the fact that the extra iso-singlet 
quarks carry zero weak isospin and thus, are not involved in SU{2)l interactions as weak 




FIG. 7: Annihilation diagram contributing to the 5+ — ^ K*'^j decay. 
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eigenstates. For example, the charge current interaction term 

5 



3 

Jcc' = Y.^^nn\^'d\W; + H.C. , (17) 



transforms to 



where 



Jcc' = E I^u'lV-^^^d'.'w; + H.C. , (18) 



y.,^J2{Air\Aif , (19) 



1=1 



when expressed in terms of mass eigenstates via eq. (fTH|l . ^ is the 4x4 generahzation of 



the Cabibbo-Kobavashi-MaskawafCKM)|l2|| quark mixing matrix. The fact that the fourth 
generation is iso-singlet, i.e. i = 1,2, 3, leads to non-unitarity of the mixing matrix V as 
demonstrated in the following: 

(5=1 5=1 ij=l 

i,j=l 5=1 

= J2 i^D^^f = - (AfyAf^ . (20) 

i=l 

In obtaining the last line in (j^Uj) . the unitarity of the tranformation matrix A^ has been 
utilized. In the same way, one can show 

Equations (j2(Jj) and (j2ip . together with the unitarity of A^"^, indicate that the quark mixing 
matrix V can not be unitary. This property of the VQM has interesting consequences in the 
neutral current sector where non-vanishing tree level FCNC, proportional to the deviation 
of the quark mixing matrix (eq. ()19p ) from unitarity, are generated. To demonstrate this 
explicitly, let us examine the neutral current which is coupled to boson 

Jnc' = (C E iLl'gl - Q, sin^ 9^ i^^qlYql + qM)] , (22) 

where Qg is the electric charge of the quark q. The first term in ()22|1 is proportional to 
the third component of the isospin, which has the value +1/2 or —1/2 for the up- or down- 
type quarks, respectively. Consequently, the iso-singlet quarks, which have zero isospin, do 
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not contribute to this term. As a result, under the transformations ()16p. the neutral current 
can be expressed in terms of the mass eigenstates as follows: 

3 4 

cos 9,, 



i=l a,P=l 
4 4 

' ,q5l3 



5=1 a,/3=l 

E (/^f/""g'L7'^gi" - sin^ 0^6-^q"^Yq") , (23) 
cosb',,, V / 



a,l3=l 



where 



1=1 



(\/t\/)"'^ , q = down — type 
{VVT^ ,q = up- type 



(24) 



We observe that the non-unitarity of the mixing matrix V in the VQM leads to the tree 
level FCNC in the Z sector. This in turn results in additional contributions to the isospin 
symmetry breaking. In fact, in the following section, we show that the isospin asymmetry in 
B K*'~f transitions offers an excellent physical observable for constraining the parameters 
of the VQM. As is shown in our result, the advantage here is that this asymmetry, unlike 
the branching ratio for B K*'~f for example, is sensitive to only one model parameter, 
namely the non-unitarity parameter U^'^. Therefore, isospin asymmetry can provide a good 
constraint on the size of the FCNC in the context of VQM irrespective of the masses of the 
additional quarks. 

IV. ISOSPIN ASYMMETRY IN B ^ k*j WITHIN THE VQM 

The non-vanishing FCNC at the tree level leads to an additional contributing Feynmann 



diagram of fig. |H1 The amplitude for bq — > sq via exchange in VQM can be written as 11| : 



2 cos(^) 



[(/^-Q,sin2^)g7^(l-75)g-QgSin2 0g7^(l + 75)g] , (25) 
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where U^'' = (V^^V^)** is a measure of the non-unitarity of the extended quark mixing matrix 
as derived in eq. |211 One can then write in terms of the effective operators O3 and 
O5 and therefore the presence of the extra vector quarks is summarized in some additional 
terms in the Wilson coefficients C3 and C5 to the leading order in the strong coupling a^. 

VtbVts 

_ U'^ r 1/2 - 2/3 sin^ 6 = 0.35 ...q = up 



VtbVtl {-1/2 + 1/3 8111^6 = -0A2...q = down ' 

^VQM _ U'' ^ 2._i^/-2/3sin2^ = -0.15...g = «p 

~ VttV:r' Vt,Vtl \l/3smH = 0.08... q = down ' ^^^^ 

With the upper bound {U'^^l < 10"'^ coming from the rare B decays [l^, the additional 
contribution due to the tree level FCNC could be comparable to the SM value of these 
coefficients at = nih, i.e. C3 = 0.014 and C5 = —0.041. 

Here an explanation is in order. Strictly speaking, one should include these extra terms, 
which are proportional to the electric charge of the light quark, in the electroweak penguin 
operators O7...io[13]. However, since, as far as isospin symmetry breaking to the leading 
order of is concerned, one can ignore these operators within SM, we prefer to write the 
additional VQM-generated contributions in terms of the dominant QCD penguin operators. 
In any case, our result does not change had we followed the strict formulation of the problem. 

Following the method of Ref. j^, the nonspectator isospin symmetry breaking contribu- 
tion can be written as Ag = bqAi^ad, with q = u 01 d being the flavor of the light anti-quark 
in the B meson. Ai^ad is the leading isospin symmetry conserving spectator amplitude. As 
mentioned in the introduction, to the leading order in the strong coupling constant a<j, the 
main contribution to -B ^ K*'-) is from the electromagnetic penguin operator and the 
factorizable amplitude A/ead is proportional to the form factor Tf^^* which parameterizes 
the hadronic matrix element of this operator to the leading order in KqcD/^b- bg is the 
parameter that depends on the flavor of the spectator. In fact, this parametrization of 
the nonspectator (nonfactorizable) contribution leads to a simple expression of the isospin 
asymmetry of eq. (jH]) in terms of bg-. 

Ao- = ^{ba - bu) , (27) 
Using the expression for bg which is derived in Ref within the QCD factorization method. 
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FIG. 8: Tree level contribution to non-spectator process in i? — > K*^. Cross represents the 
alternative coupling of the emitted photon. 

TABLE I: The numerical values of the parameters in eq. H28|) . 
Ti rrifc \b Jk' fk- '^B rriK* F± af 

0.32 4.2 GeV 0.35 GeV 0.226 GeV 0.175 GeV 5.28 GeV 0.892 GeV 1.21 -0.41 

we obtain the share of vector quarks to the isospin asymmetry as follows: 



sb 



0.38^ — ^ -0.28 



(28) 



The numerical input for the parameters of eqn. (j^H|l are tabulated in Table HI which results 
in an isospin asymmetry due to the extra generation of quarks of the form: 

Ao^_«*^ = fi^l . (29) 



The significance of eq. is that it is sensitive only to one model parameter, i.e. U''^'^, 

and therefore, with more precise experimental data becoming available in the future, this 
observable could serve to impose a stringent constrain on the all important nonunitarity 
parameter of the vector quark model. 
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V. CONCLUSION 



Isospin asymmetry in B — > i^*7 decay mode could prove to be an important observable 
in the precision test of the SM and constraining various scenarios for new physics beyond 
it. Since this asymmetry is defined as the ratio of branching fractions, its sensitivity to the 
nonperturbative parameters could be minimal. In this work, we calculated the additional 
contribution to the isospin asymmetry due to an extra generation of vector quarks and 
derived an expression in terms of the model parameter U^^ within the QCD factorization 
method. A precise measurement of BR{B K*'~f) for both neutral and charged B-mesons 
should impose a stringent constraint of this important VQM model parameter 
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